New  hyperspectral  discrimination  measure  for 
spectrai  characterization 


Yingzi  Du,  member  spie 

Laboratory  for  Biometric  Signal  Processing 

United  States  Naval  Academy 

Electrical  Engineering  Department 

Annapolis,  Maryland  21402 

E-mail:  ydu@usna.edu 

Chein-i  Chang,  fellow  spie 
University  of  Maryland,  Baltimore  County 
Department  of  Computer  Science 
and  Electrical  Engineering 
Remote  Sensing  Signal  and  Image 
Processing  Laboratory 
Baltimore,  Maryland  21250 

Hsuan  Ren,  member  spie 
National  Central  University 
Department  of  Information  Engineering 
Center  for  Space  and  Remote 
Sense  Research 
Chungli,  Taiwan 

Chein-Chi  Chang 

University  of  Maryland,  Baltimore  County 
Department  of  Civil  and  Environmental 
Engineering 

Baltimore,  Maryland  21250 


Abstract.  The  spectral  angle  mapper  (SAM)  has  been  widely  used  in 
multispectral  and  hyperspectral  image  analysis  to  measure  spectral  simi¬ 
larity  between  substance  signatures  for  material  identification.  It  has 
been  shown  that  the  SAM  is  essentially  the  Euclidean  distance  when  the 
spectral  angle  is  small.  Most  recently,  a  stochastic  measure,  called  the 
spectral  information  divergence  (SID),  has  been  suggested  to  model  the 
spectrum  of  a  hyperspectral  image  pixel  as  a  probability  distribution,  so 
that  spectral  variations  among  spectral  bands  can  be  captured  more 
effectively  in  a  stochastic  manner.  This  paper  develops  a  new  hyper¬ 
spectral  spectral  discrimination  measure,  which  combines  the  SID  and 
the  SAM  into  a  mixed  measure.  More  specifically,  let  r  and  r'  denote  two 
hyperspectral  image  pixel  vectors  with  their  corresponding  spectra  speci¬ 
fied  by  s  and  s'.  Then  SAM(s,s')  measures  the  spectral  angle  between 
s  and  s'.  Similarly,  SID(s,s')  measures  the  information  divergence  be¬ 
tween  the  probability  distributions  generated  by  s  and  s'.  The  proposed 
new  measure,  referred  to  as  the  SID-SAM  mixed  measure,  can  be  imple¬ 
mented  in  two  versions,  given  by  SID(s,s')xtan(SAM(s,s'))  and 
SID(s,s')xsin(SAM(s,s')),  where  tan  and  sin  are  the  usual  trigonomet¬ 
ric  functions.  The  spectral  discriminability  of  such  a  mixed  measure  is 
greatly  enhanced  by  multiplying  the  spectral  abilities  of  the  two  mea¬ 
sures.  In  order  to  demonstrate  its  utility,  a  comparative  study  is  con¬ 
ducted  among  the  SID-SAM  mixed  measure,  the  SID,  and  the  SAM.  Our 
experimental  results  have  shown  that  the  discriminatory  ability  of  the 
(SID, SAM)  mixed  measure  can  be  a  significant  improvement  over  the 
SID  and  SAM.  ©  2004  Society  of  Photo-Optical  Instrumentation  Engineers. 
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1  Introduction 

A  remotely  sensed  image  is  usually  collected  by  a  number 
of  spectral  channels,  each  of  which  produces  its  individual 
and  separate,  but  coregistered  image.  As  a  result,  it  is  actu¬ 
ally  an  image  cube  with  each  image  pixel  represented  by  a 
column  vector  of  which  each  component  is  a  pixel  of  an 
image  acquired  by  a  particular  spectral  channel.  With  re¬ 
cent  remote  sensing  instruments  such  as  imaging  spectrom¬ 
eters,  hundreds  of  spectral  channels  can  be  used  to  detect 
many  material  substances  that  generally  cannot  be  resolved 
by  multispectral  imaging  sensors.  Such  images  are  gener¬ 
ally  referred  to  as  hyperspectral  images,  as  opposed  to  mul¬ 
tispectral  images,  which  are  acquired  by  tens  of  spectral 
channels.  With  so  many  additional  spectral  channels  used 
for  data  acquisition,  a  hyperspectral  image  pixel  vector  pro¬ 
vides  more  spectral  information  than  does  a  multispectral 
image  pixel  vector.  In  many  situations,  such  spectral  infor¬ 
mation  is  valuable  or  even  crucial  in  data  analysis. 

In  order  to  capture  and  characterize  the  spectral  proper¬ 
ties  provided  in  a  single  pixel  vector  by  hundreds  of  spec¬ 
tral  channels,  a  new  stochastic  measure,  called  the  spectral 
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information  measure  (SIM),  was  introduced  in  Refs.  1-4.  It 
models  the  spectrum  of  a  pixel  vector  as  a  probability  dis¬ 
tribution  so  that  the  spectral  properties  of  the  pixel  vector 
can  be  further  characterized  by  its  statistical  moments  of 
any  order,  such  as  the  mean,  variance,  skewness,  and  kur- 
tosis.  More  importantly,  by  virtue  of  the  SIM,  many  con¬ 
cepts  in  information  theory  are  readily  applied  to  spectral 
characterization.  For  example,  the  self-information  derived 
from  the  SIM  can  be  used  to  describe  the  information  pro¬ 
vided  by  a  particular  spectral  channel  within  a  pixel  vector. 
Using  such  self-information,  a  discrimination  measure, 
called  the  spectral  information  divergence  (SID),  can  be 
derived  and  used  to  measure  the  spectral  similarity  between 
two  pixel  vectors. 

In  the  remote  sensing  community,  the  spectral  angle 
mapper  (SAM)  has  been  widely  used  as  a  spectral  similar¬ 
ity  measure  for  material  identification.  It  calculates  the 
angle  between  two  spectra  and  uses  it  as  a  measure  of  dis¬ 
crimination.  Another  popular  spectral  similarity  measure. 
Euclidean  distance  (ED),  has  also  been  used  to  calculate  the 
distance  between  two  spectra  as  a  spectral  similarity  mea¬ 
sure.  It  is  shown  in  Refs.  1  and  2  that  when  the  angle  is 
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small,  the  SAM  and  ED  yield  very  close  results  and  they 
are  essentially  equivalent  in  terms  of  spectral  discrimina¬ 
tion.  In  contrast  with  the  SID,  the  SAM  and  the  ED  are 
deterministic  measures  and  consider  a  spectrum  of  a  pixel 
vector  as  a  vector  rather  than  a  probability  distribution  as 
modeled  by  the  SIM. 

This  paper  develops  a  new  measure  that  combines  the 
SID  and  SAM  to  form  a  mixed  measure,  called  the  SID- 
SAM  mixed  measure,  which  takes  advantages  of  strengths 
of  both  measures.  It  can  be  implemented  in  two  versions, 
SID(s,s')Xtan(SAM(s,s'))  and  SID(s,s') 

X  sin(SAM(s,s')),  where  s  and  s'  are  the  spectra  of  the 
two  pixel  vectors  r  and  r',  and  tan  and  sin  are  the  tangent 
and  sine  trigonometric  functions  respectively.  The  SID- 
SAM  mixed  measure  takes  advantage  of  the  strengths  of 
both  the  SID  and  the  SAM  in  spectral  discriminability,  in 
the  sense  that  the  spectral  similarity  and  dissimilarity  re¬ 
sulting  from  the  mixed  measure  is  considerably  enhanced 
by  multiplying  the  spectral  abilities  of  the  two  measures. 
The  reason  for  taking  tangent  or  sine  rather  than  cosine  is 
to  calculate  the  perpendicular  distance  between  two  vectors 
instead  of  the  projection  of  one  vector  along  the  other  vec¬ 
tor.  As  a  result,  experimental  results  demonstrate  that  the 
simple  SID- SAM  mixed  measure  can  be  a  significant  im¬ 
provement  in  discriminatory  ability  over  the  SID  and  the 
SAM. 

The  remainder  of  this  paper  is  organized  as  follows. 
Section  2  briefly  reviews  the  concept  of  spectral  informa¬ 
tion  measure  developed  in  Refs.  1  and  2.  Section  3  intro¬ 
duces  three  spectral  discrimination  measures:  the  SID,  the 
SAM,  and  the  proposed  SID-SAM  mixed  measure.  Section 

4  describes  some  measures  that  can  be  used  to  evaluate  the 
effectiveness  of  a  spectral  discrimination  measure.  Section 

5  presents  experimental  results.  Section  6  concludes  with 
some  remarks. 

2  Spectral  Information  Measure 

The  spectral  information  measure  (SIM)  is  a  newly  devel¬ 
oped  stochastic  measure,^  which  considers  the  spectral 
band-to-band  variability  as  a  result  of  uncertainty  incurred 
by  randomness.  It  models  the  spectral  values  of  each  hy¬ 
perspectral  image  pixel  vector  as  a  random  variable  with 
the  probability  distribution  obtained  by  normalizing  its 
spectral  histogram  to  unity.  With  this  interpretation,  the 
SIM  is  a  measure  of  the  spectral  variability  of  a  single 
hyperspectral  image  pixel  vector  based  on  interband  corre¬ 
lation.  It  not  only  can  capture  the  randomness  of  interband 
spectral  changes  of  a  pixel  vector,  but  also  can  generate 


SpEcCroj  of  bJftdLli'Qfh:,  cteofoIe  Les-vei,  JEd  sal  ands^Ebruib 


Fig.  1  Spectra  of  five  AVIRIS  reflectances:  (a)  Spectral  signatures 
of  blackbrush,  creosote  leaves,  dry  grass,  red  soil,  and  sagebrush, 
(b)  Spectral  signature  of  sagebrush,  used  as  the  reference  signa¬ 
ture. 


high-order  statistics  of  spectral  variations.  Therefore,  the 
SIM  can  be  considered  as  a  single-pixel-based  information- 
theoretic  measure. 

Eor  a  given  hyperspectral  pixel  vector  r 
=  (rj  ,r2,...,r^)^,  each  component  rj  represents  a  pixel  in 
the  band  image  Bj ,  which  is  acquired  by  a  certain  wave¬ 
length  (x)j  in  a  specific  spectral  range.  Let  s 
=  be  the  corresponding  spectral  signature 

(i.e.,  spectrum)  of  r,  where  Sj  represents  the  spectral  signa¬ 
ture  of  r,  in  the  form  of  either  radiance  or  reflectance  val- 
ues.  Suppose  that  {o)j}j^^  is  a  set  of  L  wavelengths,  each 
of  which  corresponds  to  a  spectral  band  channel.  Then  r 
can  be  modeled  as  a  random  variable  by  defining  an  appro- 


Table  1  Discrimination  values  produced  by  SAM  and  SID. 


SAM 

Blackbrush 

Creosote  leaves 

Dry  grass 

Red 

soil 

Sagebrush 

SID 

Blackbrush 

0 

0.1767 

0.2575 

0.4058 

0.0681 

Creosote  leaves 

0.0497 

0 

0.4213 

0.5714 

0.1289 

Dry  grass 

0.0766 

0.2298 

0 

0.2179 

0.2968 

Red  soil 

0.1861 

0.4154 

0.0640 

0 

0.4515 

Sagebrush 

0.0063 

0.0303 

0.0973 

0.2340 

0 
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Table  2  Discrimination  values  produced  by  SID(SIN)  and  SID(TAN). 


SID(SIN) 

Blackbrush 

Creosote  leaves  Dry  grass 

Red  soil 

Sagebrush 

SID(TAN)  Blackbrush 

0 

0.0087 

0.0195 

0.0735 

0.0004 

Creosote  leaves 

0.0089 

0 

0.094 

0.2247 

0.0039 

Dry  grass 

0.0202 

0.103 

0 

0.0138 

0.0285 

Red  soil 

0.0800 

0.2671 

0.0142 

0 

0.1021 

Sagebrush 

0.0004 

0.004 

0.0298 

0.1135 

0 

priate  probability  space  associated  with  it,  where 

n  is  a  sample  space,  2  is  an  event  space,  and  P  is  a  prob¬ 
ability  measure.  In  our  case,  H  =  is  the 

sample  space,  2  is  the  power  set  of  II  (i.e.,  the  set  of  all 
subsets  of  H),  and  r((x)i)  =  Sj . 

In  order  to  define  a  legitimate  probability  measure  P  for 
r,  we  first  assume  that  all  components  Sj  associated  with  r 
are  nonnegative.  This  is  generally  a  valid  assumption  due  to 
the  nature  of  radiance  or  reflectance.  With  this  assumption, 
we  can  normalize  Sj  to  the  range  of  [0,  1]  as  follows. 


(1) 


Using  Eq.  (1),  we  can  further  define  a  probability  measure 
P  for  r  by 


Pi{^j})=Pj-  (2) 

The  probability  vector  P=(Pi,P2,---,Pl)^  is  the  probabil- 
ity  mass  function  of  the  probability  measure  P  and  is  the 
desired  probability  distribution  of  the  pixel  vector  r.  By 
means  of  this  probability  interpretation,  any  pixel  vector  r 
=  can  be  viewed  as  a  single  information 

source  with  its  statistics  governed  by  V=(Pi,P2^---^Pl)^ 
via  Eqs.  (1)  and  (2).  As  a  result,  p=  (P\,P2^--^PlV  can  be 
used  to  describe  the  spectral  variability  of  a  pixel  vector 
and  its  statistics  of  any  order  (mean,  variance,  skewness, 
kurtosis,  etc.).  Eor  instance,  we  can  define  its  statistics  of 
different  orders,  such  as  the  mean  =  ,  vari¬ 

ance  a^(r)  =  ^f^iPi[si— third  central  moment 
K^(r)  =  'Zf^iPi[si— ,  fourth  central  moment  Ac^(r) 
=  etc. 

Erom  information  theory^  we  can  further  use  p 
=  (Pi,P2,... ,Pl) ^  to  define  the  self-information  provided 
by  a  particular  band,  say  band  j,  by 


Ij{r)=-logpj,  (3) 

which  describes  how  much  information  is  yielded  by  the 
band  image  Bj  .  The  entropy  of  a  hyperspectral  image  pixel 
vector  r,  denoted  by  H(r),  is  actually  the  mean  of  the 
self-information  over  all  bands  and  can  be  calculated  by 


L  L 

H{r)  =  ^  lj{r)pj  =  -  2  Pj  log Pj ,  (4) 

7=1  7=1 

which  measures  the  uncertainty  resulting  from  the  pixel 
vector  r. 

3  Spectral  Discrimination  Measures 

In  this  section,  three  spectral  discrimination  measures  are 
presented  and  used  to  measure  the  similarity  between  any 
two  pixel  vectors.  Assume  pixel  vectors  r  and  r'  with  their 
respective  spectral  signatures  given  by  s= 
and  s'  =  (5j  ,s'2,...,s[Y. 

3.1  Spectral  Information  Divergence^ 

Let  p=(pi,P2,...,PlV  and  q=iqi,q2,...,qLV  be  the 
two  probability  mass  functions  generated  by  s 
=  (si,S2,...,Sl)^  and  s'  =  (s[  ,^'2  the  spectral  sig¬ 

natures  of  r  and  r',  respectively.  So  the  self-information 
provided  by  r'  for  band  j  is  defined  by  Eq.  (3)  and  given 
by 

Ij{r')=-logqj.  (5) 

Using  Eqs.  (3)  and  (5),  the  discrepancy  in  the  self¬ 
information  of  the  band  image  Bj  in  r  relative  to  the  self¬ 
information  of  Bj  in  r',  denoted  by  Dj(r\\r'),  can  be  de¬ 
fined  as 


£)/  r||  r' )  =  //  r)  -  //  r' )  =  log{pj  Iqj).  (6) 

Averaging  T)^(r||r')  in  Eq.  (6)  over  all  the  band  images 
with  respect  to  r  results  in 


Fig.  2  Graphical  plots  of  Tables  1  and  2. 
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Table  3  RSDPW  values  of  SAM  and  SID  using  sagebrush  as  the  reference  signature. 


SAM 


Blackbrush 

Creosote  leaves 

Dry  grass 

Red  soil 

SID 

Blackbrush 

0 

1 .892  805 

4.358  297 

6.629  956 

Creosote  leaves 

4.809  524 

0 

2.302  56 

3.502  715 

Dry  grass 

15.444  44 

3.211  221 

0 

1.521226 

Red  soil 

37.142  86 

7.722  772 

2.404  933 

0 

Table  4  RSDPW  values  of  SID(SIN)  and  SID(TAN)  using  red  soil  as  the  reference  signature. 


SID(SIN) 

Blackbrush  Creosote  leaves  Dry  grass  Red  soil 


SID(TAN)  Blackbrush  0 

Creosote  leaves  9. 1 72  358 

Dry  grass  68.164  23 

Red  soil  257.6992 


9.162  602  67.494  31  243.15  28 

0  7.366  282  26.537  53 

7.431  484  0  3.602  568 

28.0952  3.780  563  0 


Fig.  3  Graphical  plots  of  Tables  4  and  5. 


Table  5  RSDPB  and  self-information  values  generated  by  SAM,  SID,  SID(SIN),  and  SID(TAN). 


Measure 

Quantity 

Value 

Blackbrush 

Creosote  leaves 

Dry  grass 

Red  soil 

Sagebrush 

Entropy 

SAM 

RSDPB 

0.2266 

0.3446 

0.1074 

0.0624 

0.2589 

Self-information 

0.485  33 

0.5297 

0.3457 

0.2497 

0.5047 

2.1151 

SID 

RSDPB 

0.1879 

0.4953 

0.0549 

0.0132 

0.2486 

Self-information 

0.453  21 

0.5021 

0.2299 

0.0824 

0.4992 

1 .7668 

SID(SIN) 

RSDPB 

0.1519 

0.5953 

0.0213 

0.003 

0.2284 

Self-information 

0.412  99 

0.4455 

0.1183 

0.0251 

0.4866 

1 .4885 

SID(TAN) 

RSDPB 

0.1452 

0.6106 

0.0196 

0.0027 

0.2218 

Self-information 

0.4042 

0.434  56 

0.1112 

0.0230 

0.4819 

1 .4549 
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Fig.  4  Graphical  plots  of  self-information  values  and  entropies  in 
Table  5. 


Fig.  5  Graphical  plots  of  RSDPB  values  in  Table  5. 


L  L 

£»(r||r')  =  2  Dj{r\\r')pj='^  pjlogipj/qj),  (7) 

7=1  7=1 

where  Z)(r||r')  is  the  average  discrepancy  in  the  self¬ 
information  of  r'  relative  to  that  of  r.  In  context  of  infor¬ 
mation  theory,  Z)(r||r')  in  Eq.  (7)  is  called  the  relative  en¬ 
tropy  of  r'  with  respect  to  r;  it  is  also  known  as  the 
Kullback-Leibler  information  measure,  directed  diver¬ 
gence,  or  cross  entropy.^  Similarly,  we  can  also  define  the 
average  discrepancy  in  the  self-information  of  r  relative  to 
the  self-information  of  r'  by 


3.2  Spectral  Angle  Mapper  ^ 

The  SAM  measures  spectral  similarity  by  finding  the  angle 
between  the  spectral  signatures  s  and  s'  of  two  pixel  vec¬ 
tors,  r  and  r' : 


SAM(s,s')  =  cos  ^ 


(10) 


where  {s,s')  =  ^f=iSis'i  ,  =  and  ||s' 


£)(r'||r)=2  Dj{r'\\r)qj='2/  qjlogiqjlpj).  (8) 

7=1  7=1 

Adding  Eqs.  (7)  and  (8)  yields  the  spectral  information 
divergence  (SID),  defined  by 

SID(r,r')  =  D(r||r')  +  D(r'||r),  (9) 


3.3  SID-SAM  Mixed  Measure 

A  new  hyperspectral  measure  is  introduced  by  combing  the 
SID  and  SAM  into  a  new  measure,  referred  to  as  the  SID- 
SAM  mixed  measure,  which  can  be  implemented  in  two 
versions.  The  first  version  is  obtained  by  multiplying  the 
SID  by  the  tangent  of  the  SAM  between  two  spectral  sig¬ 
natures  s  and  s'  and  is  given  by 


which  can  be  used  to  measure  the  discrepancy  between  two 
pixel  vectors  r  and  r'  in  terms  of  their  corresponding  prob¬ 
ability  mass  functions  p  and  q.  It  should  be  noted  that  while 
SID(r,r')  is  symmetric,  D(r||r')  is  not.  That  is, 
SID(r,r')  =  SID(r',r)  and  D(r||r')  AD(r'||r). 


SID(TAN)  =  SID(s,s')Xtan(SAM(s,s')).  (11) 

The  second  version  replaces  the  tangent  function  in  Eq. 
(11)  with  the  sine  function  and  yields  the  following  new 
measure: 


(a)  (b) 


Fig.  6  (a)  A  HYDICE  panel  scene  that  contains  1 5  panels;  (b)  ground  truth  map  of  the  spatial  locations 
of  the  15  panels 
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Fig.  7  Spectra  of  P1 ,  P2,  P3,  P4,  and  P5. 


SID(SIN)  =  SID(s,s' )  X  sin(SAM(s,s' )),  (12) 

which  is  the  product  of  the  SID  and  the  sine  of  the  SAM 
between  two  spectral  signatures  s  and  s'. 

By  taking  the  product  of  two  measures,  the  SAM  and  the 
SID,  the  spectral  discriminability  of  the  new  SID-SAM 
mixed  measure  is  increased  considerably  because  it  makes 
two  similar  spectral  signatures  even  more  similar  and  two 
dissimilar  spectral  signatures  more  distinct.  However,  if  the 
cosine  function  is  used  to  replace  the  sine  function  in  Eq. 
(12),  the  spectral  discriminability  will  be  significantly  re¬ 
duced,  because  the  cosine  calculates  the  projection  of  one 
spectral  signature  along  the  other  one  instead  of  the  projec¬ 
tion  of  one  spectral  signature  orthogonal  to  the  other. 

4  Measures  of  Spectral  Discriminability 

In  this  section,  we  briefly  review  three  measures  proposed 
in  Refs.  1  and  2  that  can  be  used  to  evaluate  the  effective¬ 
ness  of  a  discrimination  measure  in  terms  of  the  spectral 
discriminatory  probability  and  power. 

4.1  Relative  Spectral  Discriminatory  Probability 
(RSDPB) 

Let  be  K  spectral  signatures  specified  by  a  given 

set  A,  which  can  be  considered  as  either  a  database  or  a 
spectral  library.  Also  assume  that  t  is  any  specific  target 
spectral  signature  to  be  identified  via  A.  We  define  the  spec¬ 
tral  discriminatory  probabilities  of  all  s  in  A  with  respect 
to  t  by 


Table  7  Discrimination  values  produced  by  SID(SIN)  and  SID(TAN). 


SID(SIN) 

PI 

P2 

P3 

P4 

P5 

SID(TAN) 

PI 

— 

0.0002 

0.0006 

0.0027 

0.0039 

P2 

0.0002 

— 

0.0001 

0.0057 

0.0076 

P3 

0.0006 

0.0001 

— 

0.0078 

0.0097 

P4 

0.0027 

0.0057 

0.0079 

— 

0.0001 

P5 

0.0039 

0.0077 

0.0098 

0.0001 

— 

pT,A(k)  = 

for  k  = 

(13) 

^i  =  i 

m{t,Si) 

where  S^im(t,s,)  is  a  normalization  constant  determined 
by  t  and  A.  The  resulting  probability  vector 

called  relative  spectral 
discriminatory  probability  (RSDPB)  of  A  with  respect  to  t, 
or  the  spectral  discriminatory  probability  vector  of  A  rela¬ 
tive  to  t. 

Using  Eq.  (13)  we  can  identify  t  by  selecting  the  mem¬ 
ber  of  A  with  the  smallest  relative  spectral  discriminability 
probability.  If  there  is  a  tie,  any  tied  member  can  be  used 
for  t. 


4.2  Relative  Spectral  Discriminatory  Entropy 
(RSDE) 

Since  =  given  by  Eq. 

(13)  is  the  relative  spectral  discriminatory  probability  vec¬ 
tor  of  t  using  a  selective  set  of  spectral  signatures,  A 
=  {Sk}f=i  ^  we  can  further  define  the  relative  spectral  dis¬ 
criminatory  entropy  (RSDE)  of  the  spectral  signature  t  with 
respect  to  the  set  A,  denoted  by  i7R3j)E(t,A),  as  follows: 

K 

^RSDE(f’^)  “  “  2 

k=  1 

Equation  (14)  provides  a  measure  of  the  uncertainty  of 
identifying  t  by  using  A  =  {sjf=i .  A  larger 
may  have  a  smaller  chance  of  identifying  t. 

4.3  Relative  Spectral  Discriminatory  Power 
(RSDPW) 

Assume  that  m(  ,  )  is  any  given  hyperspectral  measure 
such  as  the  SID,  SAM,  or  SID-SAM  mixed  measure.  Let  d 


Table  6  Discrimination  values  produced  by  SAM  and  SID. 


SAM 

PI 

P2 

P3 

P4 

P5 

SID 

PI 

— 

0.0435 

0.0673 

0.1144 

0.1240 

P2 

0.0039 

— 

0.0430 

0.1479 

0.1567 

P3 

0.0086 

0.0033 

— 

0.1652 

0.1710 

P4 

0.0233 

0.0385 

0.0476 

— 

0.0248 

P5 

0.0313 

0.0485 

0.0570 

0.0025 

— 

Fig.  8  Graphical  plots  of  Tables  6  and  7. 
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Table  8  RSDPB  and  self-information  values  generated  by  SAM,  SID,  SID(SIN),  and  SID(TAN)  using 
P23  ■ 


Value 


Measure 

Quantity 

PI 

P2 

P3 

P4 

P5 

Entropy 

SAM 

RSDPB 

Self-information 

0.1309 

0.3840 

0.0670 

0.2613 

0.1284 

0.3802 

0.3259 

0.5271 

0.3477 

0.5299 

2.0826 

SID 

RSDPB 

Self-information 

0.0705 

0.2697 

0.0200 

0.1129 

0.0536 

0.2263 

0.3786 

0.5305 

0.4773 

0.5093 

1 .6487 

SID(SIN) 

RSDPB 

Self-information 

0.0302 

0.1525 

0.0044 

0.0344 

0.0225 

0.1232 

0.4022 

0.5285 

0.5407 

0.4797 

1.3182 

SID(TAN) 

RSDPB 

Self-information 

0.0298 

0.1510 

0.0043 

0.0338 

0.0223 

0.1224 

0.4020 

0.5285 

0.5415 

0.4792 

1.3149 

be  the  spectral  signature  of  a  reference  pixel  vector.  Sup¬ 
pose  that  s  and  s'  are  the  spectral  signatures  of  any  pair  of 
pixel  vectors.  The  RSDPW  of  m(  ,  ),  denoted  by 
RSDPW'''(s,s';d),  is  then  defined  by 

f  m(s,d)  m(s',d)] 

RSDPW''^(s,s';d)  =  max— (15) 
[m(s  ,d)  m(s,d)  J 

More  precisely,  RSDPW"'(s,s' ;d)  calculates  two  ratios — 
the  ratio  of  m(s,s';d)  to  m(s',s;d)  and  the  ratio  of 
m(s',s;d)  to  m(s,s  ';d)  — and  selects  the  greater  of  them  as 
the  discriminatory  power  of  m(  ,  ).  The  function 
RSDPW'” (s, s' ;d)  defined  by  Eq.  (15)  provides  a  quantita¬ 
tive  index  of  the  spectral  discrimination  ability  of  a  specific 
hyperspectral  measure  m(  ,  )  between  two  spectral  signa¬ 
tures  s  and  s'  relative  to  d.  Obviously,  the  higher 
RSDPW'”(s,s'  ;d)  is,  the  better  the  discriminatory  power  of 
m{  ,  )  is.  In  addition,  RSDPW'”(s,s' ;d)  is  symmetric  and 
bounded  below  by  1,  i.e.,  RSDPW'”(s,s' ;d)^  1,  with 
equality  if  and  only  if  s=s'. 

5  Experiments 

Two  data  sets  were  used  for  experiments.  The  first  was  the 
AVIRIS  (Airborne  Visible/Infrared  Imaging  Spectrometer) 
reflectance  data  shown  in  Fig.  1(a),  consisting  of  five  field 
reflectance  spectra  (blackbrush,  creosote  leaves,  dry  grass, 
red  soil,  and  sagebrush)  with  a  spectral  coverage  from  0.4 
to  2.5  /mm.  Only  158  bands  were  left  after  the  water  bands 
were  removed.  Using  these  five  signatures,  we  calculated 


PI  P2  P3  P*  P5  Entropy 

Fig.  9  Graphical  plots  of  RSDPB  values  in  Table  8. 


the  discrimination  values  for  four  spectral  measures:  the 
SAM,  the  SID,  the  SID(TAN),  and  the  SID(SIN).  Table  1 
tabulates  the  results  for  the  SAM  and  the  SID,  where  the 
values  in  the  upper  triangular  matrix  were  generated  by  the 
SAM  and  the  values  in  the  lower  triangular  matrix  were 
produced  by  the  SID.  Similarly,  the  values  generated  by  the 
SID(SIN)  and  the  SID(TAN)  are  tabulated  in  the  upper  and 
lower  triangular  matrices  of  Table  2,  respectively.  It  is  dif¬ 
ficult  to  see  which  measure  is  more  effective,  merely  from 
a  comparison  between  Table  1  and  Table  2.  Thus,  Fig.  2 
also  plots  the  values  in  Tables  1  and  2  for  visual  compari¬ 
son.  However,  despite  the  fact  that  the  SID(TAN)  and  the 
SID(SIN)  yielded  smaller  values  than  those  produced  by 
the  SAM  and  the  SID,  we  cannot  conclude  that  a  larger 
spectral  similarity  value  means  better  discriminatory  power. 

In  order  to  remedy  this  diffuculty,  we  calculate  the  rela¬ 
tive  spectral  discriminatory  power  defined  by  Eq.  (15)  us¬ 
ing  sagebrush  as  the  reference  signature  d  as  shown  in  Fig. 
1(b).  We  selected  sagebrush  for  this  purpose  because  sage¬ 
brush  is  very  close  to  both  blackbrush  and  creosote  leaves, 
as  shown  in  Tables  1  and  2,  and  provides  a  good  case  to 
evaluate  the  effectiveness  of  a  spectral  measure.  Tables  3 
and  4  tabulate  their  results.  As  we  can  see  from  these  two 
tables,  the  proposed  SID(TAN)  and  SID(SIN)  performed 
significantly  better  than  the  SID  and  the  SAM.  Figure  3 
plots  their  RSDPW  values,  which  also  clearly  indicate  the 
superior  performance  of  the  SID(TAN)  and  the  SID(SIN)  to 
the  SID  and  the  SAM.  For  example,  the  spectral  signatures 
of  sagebrush,  blackbrush,  and  creosote  leaves  are  very 


Fig.  10  Graphical  plots  of  self-information  values  and  entropies  in 
Table  8. 
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Table  9  RSDPW  values  of  SAM  and  SID  using  P2  as  the  reference 
signature. 


SAM 

PI 

P3 

P4 

P5 

SID 

PI 

— 

1.0116 

3.4000 

3.6023 

P3 

1.1818 

— 

3.4395 

3.6442 

P4 

9.8718 

11.6667 

— 

1 .0595 

P5 

12.4103 

14.6667 

1.2571 

— 

Table  10  RSDPW  values  of  SID(SIN)  and  SID(TAN)  using  P2  as 
the  reference  signature. 


SID(SIN) 

PI 

P3 

P4 

P5 

SID(TAN) 

PI 

— 

2.0000 

28.5000 

38.0000 

P3 

2.0000 

— 

57.0000 

76.0000 

P4 

28.5000 

57.0000 

— 

1 .3333 

P5 

38.5000 

77.0000 

1 .3509 

— 

close  according  to  Fig.  1.  As  expected,  the  spectral  dis¬ 
crimination  among  these  three  is  difficult.  Apparently, 
based  on  our  experiments,  the  SAM  was  the  worst,  with 
discriminatory  power  1.892  805.  The  SID  improved  on  the 
discriminatory  power  of  the  SAM  by  2.5  times;  for  the 
SID-SAM  mixed  measure  the  discriminatory  power  was 
almost  5  times  that  of  the  SAM  and  twice  that  of  the  SID. 
It  is  worth  noting  that  the  RSDPW  values  between  dry 
grass  and  red  soil  are  very  small  among  all  the  four  mea¬ 
sures.  This  is  mainly  because  the  sagebrush,  which  was 
selected  as  the  reference  signature,  has  a  very  distinct  spec¬ 
tral  signature  from  those  of  dry  grass  and  red  soil.  As  a 
result,  the  RSDPW  values  yielded  by  the  four  measures 
were  relatively  small,  which  implies  that  the  discriminatory 
power  between  dry  grass  and  sagebrush  was  nearly  the 
same  as  that  between  red  soil  and  sagebrush.  These  results 
also  suggested  that  the  dry  grass  and  red  soil  have  very 
distinct  signatures  from  the  signature  of  sagebrush. 

In  order  to  further  evaluate  the  ability  of  the  mixed  mea¬ 
sures  in  discriminatory  power  against  a  reference  signature, 
we  used  the  same  target  signature  used  in  Ref.  1,  which  is 
a  mixed  signature  composed  of  0.1055  blackbrush,  0.0292 
creosote  leaves,  0.0272  dry  grass,  0.7588  red  soil,  and 
0.0974  sagebrush.  Table  5  tabulates  the  RSDPB  and  self¬ 
information  values  with  their  entropies  for  the  signatures; 
the  proposed  SID-SAM  mixed  measures,  SID(SIN)  and 
SID(TAN),  yielded  the  two  smallest  entropies.  Plots  of 
comparisons  among  these  four  measures  are  also  shown  in 
Figs.  4  and  5,  where  the  proposed  SID(SIN)  and  SID(TAN) 
clearly  have  better  spectral  discriminability  in  identifying 
the  target  signature  t  as  red  soil,  with  smaller  spectral  dis¬ 
criminatory  probabilities  than  those  produced  by  the  SID 
and  the  SAM. 

The  data  in  the  second  set  to  be  used  in  the  experiments 
are  real  HYDICE  (Hyperspectral  Digital  Image  Collection 
Experiment)  data.  A  HYDICE  image  scene  is  shown  in  Eig. 
6(a)  (band  80)  with  the  ground  truth  map  in  Eig.  6(b).  The 
spatial  positions  of  these  15  panels  are  precisely  located, 
with  black  (B)  pixels  indicating  the  panel  center  pixels  of 
all  the  15  panels,  and  white  (W)  pixels  being  panel  pixels 
mixed  with  background  pixels.  In  this  case,  a  W  pixel  can 
be  considered  as  a  mixed  pixel.  The  scene  has  size  64 
X64  pixels;  the  low-signal,  high-noise  bands  (bands  1  to  3 
and  bands  202  to  210),  and  water  vapor  absorption  bands 
(bands  101  to  112  and  bands  137  to  153)  have  been  re¬ 
moved.  The  15  panels  are  arranged  in  a  5  X  3  matrix.  Each 
element  in  this  matrix  is  a  square  panel  and  denoted  by  pij 
with  row  indexed  by  /=!,..., 5  and  column  indexed  by  j 
=  1,2,3.  Eor  each  row  /=!,..., 5,  the  three  panels 


Pn,Pi2^Pi3  were  painted  with  the  same  material  but  have 
three  different  sizes.  Eor  each  column  7  =  1,2,3,  the  five 
panels  P\j  ,P2j  ^P3j  ^Psj  have  the  same  size  but  were 
painted  by  five  different  materials.  The  sizes  of  the  panels 
in  the  first,  second,  and  third  columns  are  3mX3m,  2m 
X2m,  and  ImX  Im,  respectively.  The  15  panels  have  five 
different  materials  and  three  different  sizes.  Eigure  7  plots 
the  five  panel  spectral  signatures  obtained  from 

Eig.  6(b),  where  the  /’th  panel  signature,  denoted  by  P/, 
was  generated  by  averaging  all  B  pixels  in  row  I  and  was 
used  to  represent  target  knowledge  of  the  panels  in  row  /. 
The  1.5-m  spatial  resolution  of  the  image  scene  suggests 
that  except  for  P2\^P3\^P4\^P5\^  which  are  two-pixel  pan¬ 
els,  all  the  rest  of  the  11  panels  are  single-pixel  panels. 

As  with  the  AVIRIS  data  studied  previously,  we  con¬ 
ducted  similar  experiments  using  these  five  panel  signa¬ 
tures.  Tables  6  and  7  tabulate  spectral  discrimination  values 
of  the  SAM,  the  SID,  the  SID(SIN),  and  the  SID(TAN);  the 
values  in  the  upper  triangular  matrices  were  generated  by 
the  SAM  and  the  SID(SIN),  respectively,  while  the  values 
in  the  lower  triangular  matrices  were  produced  by  the  SID 
and  the  SID(TAN),  respectively.  As  expected,  the  values  of 
the  proposed  SID(SIN)  and  SID(TAN)  were  also  smaller 
than  those  produced  by  the  SAM  and  the  SID.  In  analogy 
with  Pig.  2,  Pig.  8  plots  the  values  in  Tables  6  and  7  for 
graphical  comparison.  Once  again,  we  also  calculated  the 
RSDPB  and  self-information  values  among  the  five  panel 
signatures  and  their  entropies.  The  results  are  tabulated  in 
Table  8,  with  their  corresponding  graphical  plots  in  Pigs.  9 
and  10.  As  shown  in  these  tables  and  Pig.  9,  the  proposed 
SID(SIN)  and  SID(TAN)  again  produced  the  smallest  en¬ 
tropies.  The  RSDPW  values  of  the  four  measures  are  tabu¬ 
lated  in  Tables  9  and  10,  and  their  corresponding  graphical 
plots  are  in  Pig.  11,  where  the  reference  signature  d  was 
chosen  to  be  the  panel  signature  P23.  As  shown,  the  RS¬ 
DPW  values  of  the  proposed  SID(SIN)  and  SID(TAN) 
yielded  better  discriminatory  powers  than  did  the  SID  and 


Fig.  11  Graphical  plots  of  RSDPW  values  in  Tables  9  and  10. 
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the  SAM.  The  magnitude  of  the  improved  power  is  propor¬ 
tional  to  the  dissimilarity  between  the  signatures  and  the 
reference  signature. 


6  Conclusions 

This  paper  presents  a  new  hyperspectral  discrimination 
measure,  the  SID-SAM  mixed  measure,  which  is  a  mixture 
of  the  widely  used  spectral  angle  mapper  (SAM)  and  a 
recently  developed  information  measure,  the  spectral  infor¬ 
mation  divergence  (SID).  Despite  the  fact  that  the  SID- 
SAM  mixed  measure  is  obtained  by  simply  multiplying  the 
SID  by  the  SAM,  the  improvement  on  the  performance  in 
spectral  discriminability  is  remarkable,  as  demonstrated  in 
experiments.  As  a  result,  two  similar  spectral  signatures  are 
made  more  similar,  while  two  dissimilar  spectral  signatures 
become  more  distinct.  This  is  because  the  measure  com¬ 
bines  the  strengths  of  SAM  and  SID. 
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